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OO : ABSTRACT 

. . Context. Galactic Cepheids in the vicinity of the Sun have a residual line-of-sight velocity, or 7-velocity, which shows 

' a systematic blueshift of about 2 km s~^compared to an axisymmetric rotation model of the Milky Way. This term is 

either related to the space motion of the star and, consequently, to the kinematic structure of our Galaxy, or it is the 
^ , result of the dynamical structure of the Cepheids' atmosphere. 

■ Aims. We aim to show that these residual 7- velocities are an intrinsic property of Cepheids. 

' Methods. We observed nine galactic Cepheids with the HARPS** spectroscope, focusing specifically on 17 spectral 

I lines. For each spectral line of each star, we computed the 7-velocity (resp. 7-asymmetry) as an average value of the 

interpolated radial velocity (resp. line asymmetry) curve. 

Results. For each Cepheid in our sample, a linear relation is found between the 7-velocities of the various spectral 
lines and their corresponding 7-asymmetries, showing that residual 7-velocities stem from the intrinsic properties of 
Cepheids. We also provide a physical reference to the stellar 7-velocity: it should be zero when the 7-asymmetry is zero. 
Following this definition, we provide very precise and physically calibrated estimates of the 7-velocities for all stars of 
. our sample [in km s"^]: -11.3±0.3 [R TrA], -3.5±0.4 [S Cru], -1.5±0.2 [Y Sgr], 9.8±0.1 [/3 Dor], 7.1±0.1 [( Gem], 

■ 24.6±0.4 [RZ Vel], 4.4±0.1 [£ Car], 25.7±0.2 [RS Pup]. Finally, we investigated several physical explanations for these 
7-asymmetries like velocity gradients or the relative motion of the line-forming region compared to the corresponding 

, mass elements. However, none of these hypotheses seems to be entirely satisfactory to explain the observations. 

(«~^ . Conclusions. To understand this very subtle 7-asymmetry efi'ect, further numerical studies are needed. Cepheids' at- 

QQ ' mosphere are strongly affected by pulsational dynamics, convective flows, nonlinear physics, and complex radiative 

^— ■») ' transport. Hence, all of these effects have to be incorporated simultaneously and consistently into the numerical models 

to reproduce the observed line profiles in detail. 

Key words. Techniques: spectroscopic - Stars: atmospheres - Stars: oscillations (including pulsations) - (Stars: vari- 
ables): Cepheids - Stars: distances 



- 1 — I 



1. Introduction ing the kinematics of the Milky Way might be closely re- 
lated to the dynamical structure of Cepheids' atmosphere. 

Cepheids are very important astrophysical objects due to Concerning the distance scale, near-infrared interferom- 

their well-known period-lummosity (PL) relation. Based on ^^^.y currently provides a new, quasi-geometrical way to de- 

this relation, a multi-decade work allowed us to determine termine the distance of galactic Cepheids up to 1 kpc (see 

the kinematic structure of the Milky Way (m particular ^ g_ Sassclov & Karovska Hool and KcrvcUa et al. ((fool ). 

Its rotation) and to reach cosmologically significant extra- ^hc basic principle of the Interferomctric Baade-Wesselink 

galactic distances (see Hubble Space Telescope Key Project, ^^^-^^^ (jbw) is to compare the Hnear and angular size 

Freedman et al. 2001). In the second paper of this se- yamation of a pulsating star in order to derive its distance 

ries, Nardetto et al. (2007) (hereafter Paper II) estabhshed through a simple division. The key point is that intcrfero- 

a clear link between the distance scale and the dynami- ^^^^^^ measurements in the continuum lead to angular di- 

cal structure of Cepheids atmosphere through a period- ameters corresponding to the photosphemc layer, while the 

projection factor (Pp) relation. Similarly, studies concern- 1^^^^^^^ g^g^^^ ^^^-^^ variation is deduced by spectroscopy, 

i.e., based on line-forming regions which form higher in the 

* Based on observations made with ESO telescopes at the atmosphere. Thus, radial velocities Vrad, which are derived 

SiUa Paranal Observatory under programme IDs 072.D-0419 and from line profiles, include the integration in two directions: 

073.D-0136 over the stellar surface through limb-darkening and over 

** High Accuracy Radial velocity Planetary Search project de- the atmospheric layers through velocity gradients. All these 

veloped by the European Southern Observatory phenomena are currently merged in one specific quantity. 
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generally considered constant with time: the projection fac- 
tor p, defined as Vpuis = pV^ad, where Vpuis is defined as the 
photospheric pulsation velocity (Nardetto et al. 2004). Vpuis 
is then integrated with time to derive the photospheric ra- 
dius variation. The precision in the distance currently ob- 
tained with the IBW method is a few percent. However, it 
remains strongly dependent on the projection factor. If a 
constant projection factor is used (generally p = 1.36 for 
all stars) to derive the PL relation, errors of 0.10 and 0.03 
on the slope and zero-point of the PL relation can be in- 
troduced. This means that distances can be overestimated 
by 10% for long-period Cepheids (Paper II). 

In Paper II, we divide the projection factor into three 
sub-concepts: (1) a geometrical effect related to the limb- 
darkening, (2) the velocity gradient within the atmosphere, 
and (3) the differential motion of the "optical" pulsating 
photosphere compared to the corresponding mass elements, 
called /o-g- Even if the Pp relation was recently confirmed 
by HST observations (Fouque et al. 2007), the /o-g is rel- 
atively uncertain and currently entirely based on hydrody- 
namical calculations. Nevertheless, a key point is that this 
quantity should be, in principle, related to the so-called 7- 
velocity term. Indeed, as already mentioned by Sabbey et 
al. (1995): 

"The observed unequal line asymmetry magnitudes dur- 
ing contraction and expansion [also observed in the first pa- 
per of this series, Nardetto et al. 2006, hereafter Paper I], 
are due to the varying depth of spectral line formation over 
a Cepheid pulsation cycle. In other words, the photospheric 
spectral lines in a Cepheid atmosphere are not associated 
with the same gas particles during cycle. Therefore, they 
are not required to, and indeed do not, comply with path 
conservation 



where the integral is over a complete Cepheid cycle. 
However path conservation is a basic assumption in the 
BW method (Gautschy 1987)." 

The state of the art of the BW methods, concerning 
the 7-velocity is the following: The 7-velocity is generally 
removed, i.e. the average value of the radial velocity curve is 
forced to zero before integrating. However, the differential 
motion between the line-forming region and the gas also 
modifies the velocity amplitude. The /o-g was introduced, 
based on hydrodynamical models in Paper II, in order to 
correct this velocity amplitude effect. The next step is to 
find a correlation between /o-g and the 7-velocity. In this 
paper, we will provide some indications, but the problem 
remains unsolved. 

The 7-velocities are also of great importance for the 
determination of the kinematic structure of the Milky 
Way. Galactic Cepheids in the solar vicinity show a 
residual line-of-sight velocity (a "K-term") in their ra- 
dial velocities which is systematically blueshifted of about 
2 km s~^ compared to an axisymmetric rotation model of 
the Milky Way (Camm 1938, 1944; Parenago 1945; Stibbs 
1956; Wielen 1974; Caldwell & Coulson 1987; Moffett & 
Barnes 1987; Wilson et al. 1991; Pont, Mayor & Burki 
1994). Wielen (1974) found no correlation of the K-term 
with any obvious parameter such as period, amplitude or 
distance, and concluded that the K-term is an intrinsic 
property of Cepheid atmospheres. On the other hand, Pont 
et al. (1994) tried to revive Camm (1944) and Parenago's 



(1945) suggestion that the K-term is due to a real effect in 
the dynamics of the Galaxy. In addition, Butler et al. (1996) 
found a 7-velocity reduced by 2 km s""'^ by introducing ve- 
locity gradients in hydrostatic stellar atmospheres models. 
The K-term problem is still a matter of debate today. 

Based on very high quality HARPS observations and 
careful methodology (Sect. [2), we will show that 7-vclocities 
are due to intrinsic properties of Cepheids (Sect. [3]). Finally, 
we discuss our results in Sect. 21 

2. HARPS observations 

Ten stars have been observed with the HARPS spectrom- 
eter {R = 120000): R Tra, S Cru, Y Sgr, f3 Dor, ( Gem, 
Y Oph, RZ Vel, £ Car, RS Pup and X Sgr. X Sgr is an atyp- 
ical Cepheid presenting several components in the spectral 
line profiles. It was studied separately by Mathias et al. 
(2006). Y Oph is not studied here in detail due to its in- 
sufficient phase coverage (see Paper II, Fig. 3). We thus 
consider 8 Cepheids in thi s pape r. 

Using Kurucz models (jl992f ) we have identified about 
150 unblended spectral lines. In Paper II, we carefully se- 
lected 17 spectral lines following two criteria: (1) the con- 
tinuum must to be perfectly defined for all pulsation phases 
of all stars, in order to avoid bias in the determination of 
the line depth, (2) the selected sample of lines has to cover 
a large range of depth. The spectral lines selected are pre- 
sented in Tab. 1 of Paper II. 

As in Paper I, we use bi-Gaussian fits to derive line 
asymmetries. We repeat here the main equations in order 
to show that there is no a priori link between the line asym- 
metry and our centroid method of the radial velocity deter- 
mination. The centroid radial velocity (RVc), or the first 
moment of the spectral line profile, has been estimated as 



lunJWdX 



(1) 



where S{X) is the observed line profile. Then, the radial ve- 
locity corresponding to the minimum pixel, full width at 
half-maximum (FWHM), and asymmetry are derived si- 
multaneously by applying a classical minimization algo- 
rithm between the observed line profile (S'(A)) and a mod- 
eled spectral line profile (/(A)). The corresponding reduced 
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with N being the number of pixels in the spectral line, v 
the number of degrees of freedom and cr(Ai) = SNR *f{Xi) 
the statistical uncertainty associated to each pixel. SNR is 
the estimate of the signal-to-noise Ratio in the continuum. 
The analytic line profile is defined by 



fa^ 1 n f 41n2(A-A„,)^ 

^(^^^^'^'^"H (™m(i + ^))^ 



and 



ff^^ ^ n / 41n2(A-A^ 

^(^^ = ^-^^"H (™m(i-A))^ 



if A > A,, 



if A < A„ 



(3) 



(4) 



with four free parameters: D, the depth of the line (dimen- 
sionless); Am, the wavelength associated to the minimum of 
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Fig. 1. Fcl 4896.4 {D ~ 8%) and Fel 6024.1 {D ~ 30%) 
spectral lines evolution of f3 Dor. Pulsation phases are given 
on the left of each profile. Wavelengths have been trans- 
lated into velocities for comparison (positive velocities cor- 
respond to a rcdshift, motion toward us). The systematic 
difference in lines asymmetries and velocities is clear. The 
solid and dashed lines at = 0.61 show that modifying arti- 
ficially the line asymmetry induces a change in the resulting 
centroid velocity (see Sect.lH). 



the line (in A); FWHM, the Full- Width at Half-Maximum 
in the line (in A); and A, the line asymmetry relative to 
the FWHM (in %). Online tables (Nardetto et al. 2007) 
present the resulting values of RVm, RVc, FWHM, D, A, 
SNR, and Xred together with the corresponding uncertain- 
ties computed from the fitting method. In the following, 
we will only use RVc and A, which are clearly independent 
given the manner of their determinations. 

We insist on the RVc definition of the radial velocity 
since it is absolutely required to allow direct comparisons 
between 7- velocities of different spectral lines from different 
Cepheids. Indeed, it is the only method which provides a 
radial velocity independent of the rotation (projected on 
the line of sight) and the natural width of the spectral line 
(Burki et al. 1982 and Paper I). 

The RVc and A quantities for 17 selected spectral lines 
for each of the 8 stars have been interpolated over the pulsa- 
tion phase using a periodic cubic spline function. The inter- 
polation is performed either directly on the observational 
points (e.g. (3 Dor) or using arbitrary pivot points (e.g. 
RZ Vel). In the latter case, a classical minimization process 
between the observations and the interpolated curve is used 
to optimize the position of the pivot points. All interpo- 
lated curves presented in this study were derived using one 
of these two methods. Finally, A^ and Vj are calculated by 
averaging the A((j)) and RVc{(j)) interpolated curves, where 
(/) is the pulsation phase. The corresponding uncertainties 
are defined as the average values of individual uncertainties 
on the observational points. 

3. A relation between 7-velocities and 
7-asymmetries 

The aim of this section is to study the 7-velocities and to 
show that they consist of two components: one related to 
the space motion of the star itself, and one (the K-term) re- 
lated to the dynamical structure of Cepheids' atmosphere. 

3.1. Methodology presented in the case of P Dor 

To introduce our methodology, we discus the exemplary 
case of (3 Dor in detail. 

First, interpolated radial velocity curves derived from 
all spectral lines were corrected by the 7-velocity found 
in the Galactic Cepheid Database of the David Dunlap 
ObservatorjQ (Fernic et al. 1995, hereafter KycCD, see 
Tab. [I]). 

Fig.[T]shows the spectral line profile of two metallic lines 
Fel 4896.4 {D ~ 8% of the continuum) and Fel 6024.1 
{D c:i 30%) as a function of the pulsation phase. We trans- 
lated wavelengths into velocities for comparison. We used 
spectral lines with different line depths for clarity, but the 
results are actually independent of the line depth. Two qual- 
itative key-points concerning Fig. [1] should be mentioned: 
(1) the Fel 4896.4 spectral line (smaller line depth) seems 
to be systematically redshifted compared to the Fel 6024.1 
spectral line (look, for instance, at the pixel minimum of 
each profile) and (2) its asymmetry is also systematically 
larger (in absolute value) from phase (j) = 0.82 to ^ = 0.23, 
while it is systematically lower (in absolute value) from 
phase (j) = 0.33 to = 0.61. This leads us to the idea of a 



http:/ /www. astro. utoronto.ca/DDO/research/cepheids/ 
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Fig. 2. RVc (a) and A (b) arc represented as a function of the pulsation phase for three spectral lines in the case of 
P Dor: Fel 4896.4 (dashed line), Fel 5373.7 (sohd line), and Fel 6024.1 (dot-dashed hne). The RVc curves include only 
the GCD 7- velocity correction. The actual measurements (crosses) are only indicated in the case of the Fe I 5373.7 line. 
The corresponding uncertainties are too small to be visible in the plot. Horizontal lines correspond to the average values 
and of the interpolated curves, respectively, (c) RVc — A plots and the corresponding {A^, Vj) average values 
(crosses) for the three different lines. Although the RVc — A plots have different shapes, the average values are aligned, 
(d) Generalization of diagram (c) for all spectral lines. The RVc — A plots are not included for clarity. The upper values 
are without any correction except the GCD 7- velocity. The origin of the plot is then used as a physical reference for all 
spectral line 7- velocities of the star (lower values). The V^A^ correlation found is a strong indication that 7- velocities 
are related to intrinsic physical properties of Cepheids' atmospheres. 



correlation between the 7- velocity and the 7-asymmetry. In 
Fig. [2^b we present the corresponding interpolated radial 
velocity and asymmetry curves for these two spectral lines. 
Another line of intermediate depth is also presented (Fcl 
5373, solid line). From the interpolated curves, we now cal- 
culate the 7-velocities and -asymmetries corresponding to 
each spectral line (horizontal lines). An anti-correlation is 
clearly seen. 

In Fig. [^t, RVc — A plots and the corresponding {A^, 
Ky) average values (crosses) for the three different lines are 
presented. Although the RVc — A plots have different shapes 
(and this is confirmed for all spectral lines) , which is the re- 
sult of the dynamical structure of the Cepheid atmosphere, 
the average values (big crosses) are aligned, confirming the 
correlation found in Fig. [2^b. 

Fig.[2]i is a generalization of diagram (c) for all spectral 
lines. The RVc~A curves are not included for clarity. Upper 
values correspond to residual 7-velocities V^{i) of the 17 



selected spectral lines i after the GCD 7- velocity correction. 
A linear fit is performed, and we find the relation 

Ky(l) - KyGCD = 00^7(1) -I- bo, 

with ao = -0.16 ± 0.02 [km s'^per %] and bo = 2.4 ± 0.1 
[km s^^]. The origin of the plot is taken as a reference: the 
7-velocity is assumed to be zero when the 7-asymmetry is 
zero. This means that all points 

{A^{i),V^{i) -V^gcb) 

are translated into 

{A^{i),V^{i) - {V^GCD + bo)) 

which allows the definition of a physically calibrated 7- 
velocity for (3 Dor 

V-y^ = KyGCD + bo = 7.4 + 2.4 = 9.8±o.i km s'^ . 
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Fig. 3. Same as Fig. [2jl for all Cepheids of our sample. The origin of the VyA^ plots are used as physical references. The 
correction applied (Ky* = T^GCD + bo), together with the slope values (ao), are indicated in Tab. [1] In the last panel an 
weighted average is done over all stars. The resulting V-yA^ linear curve is very precise (ao = 0.15 ± 0.01). 



This quantity is very precise due to the very high S /N and 
spectral resolution of HARPS data. 

In principle, the line asymmetry and the 7-asymmetry 
are supposed to be the result of the dynamical structure 
of the Cepheid atmosphere only. The Vy{i)Ay{i) (here- 
after VyA^) correlation found is then a strong indication 
that residual 7-velocities (after correction by the GCD 7- 
velocities) are related to intrinsic physical properties of 
Cepheids' atmosphere and not to a real effect in dynam- 
ics of the Galaxy (see discussion). 

3.2. VyAj linear curves of all stars of our sample 

V~f and A~f were derived for all spectral lines and for all 
stars of our sample using the same method as presented in 
the case of /3 Dor. 

Tab. [1] gives V^gcd, the slope ao of the interpolated 
curves, the correction ho applied, and our final 7- 
velocities Ky* for each star. The V^A^ plots are shown in 
Fig. [21 Linear correlation curves between and A-^ are 
found for all stars of our sample, and our 60 corrections 
range from 0.2 to 3.6 km s~^. The average value is 1.8 ±0.2 
km s~^, which is consistent with the 2 km s~^"K-term" 



found in the literature. We discuss these important obser- 
vational results in the next section. 



4. Discussion 

4.1. ^-velocities 

The linear relation between and can be easily un- 
derstood. The basic principle is demonstrated in Fig. [1] for 
the pulsation phase 4> = 0.61. For clarity, we only present 
the argument for one pulsation phase. 

The solid line is the bi-Gaussian fit of the Fel 6024.1 
spectral line. Artificially decreasing the bi-Gaussian asym- 
metry of this spectral line by the average asymmetry of 
the Fel 6024.1 line, = 7.58%, as shown in Fig. [1] with 
a dashed-line, we find that the centroid-velocity increases 
(i.e. is redshifted) by an amount of -1-1.45 km s^^. To il- 
lustrate this, we included in Fig. [1] the computed line po- 
sitions RVc by vertical lines before (solid line) and after 
(dot-dashed line) modification of the bi-Gaussian asymme- 
try. Generally, we find that forcing the average line asym- 
metry Ay to zero by uniformly changing the line asymme- 
try hy the same amount at all phases alters the derived 
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Table 1. Linear relations between the 7- velocities and the 7-asymnietries Vy = uqA^ + 60 for all stars. The subscripts 

give the la uncertainty. The reduced x^, defined as Xred = with N being the number of spectral lines and v the 
number of degrees of freedom, is also indicated. The radial velocities are initially corrected from the GCD 7-velocities, 
and then a second correction bo is applied, taking the origin of A^Vy plots as a reference: = implies — 0. The 
last column corresponds to the final physically calibrated 7-velocities for each star. 



Name 


HD 


p(a) 


fflo 


Xrcd 




correction bo 








[days] 


[km s-V%] 




[km s-i] 


[km s-i] 


[km s-i] 


R TrA 


135592 


3.38925 


— 0.13±o.o5 


3 


-13.2 


1.9±0.3 


-11.3±o.3 


S Cru 


112044 


4.68976 


— 0.16±o.o7 


2 


-7.1 


3.6±o.4 


— 3.5±o.4 


Y Sgr 


168608 


5.77338 


-0.19±o.o4 


14 


-2.5 


1.0±0.2 


-1.5±0.2 


(3 Dor 


37350 


9.84262 


— 0.16±o.o2 


30 


7.4 


2.4±o.i 


9.8±o.i 


C Gem 


52973 


10.14960 


— 0.11±0.02 


22 


6.9 


0.2±o.i 


7.1±o.i 


RZ Vel 


73502 


20.40020 


-0.14±o.o9 


1 


24.1 


0.5±o.4 


24.6±o.4 


£ Car 


84810 


35.55134 


-0.15±o.o3 


15 


3.6 


0.8±o.i 


4.4±o.i 


RS Pup 


68860 


41.51500 


— 0.13±o.o5 


4 


22.1 


3.6±o.2 


25.7±o.2 



" The corresponding Julian dates (To) can be found in Paper II. 



centroid velocities RVc in such a way that the resulting 7- 
velocity becomes zero. For instance, the point {A-^ = 7.58%, 
Vy = 1.45 km s^^) of the V^A-^ linear curve of /3 Dor (Fig. 
[3|) translates into {A^ = 0%, Vy ~ km s^^). 

This result is obtained universally, regardless which 
spectral line or star is considered. Basically, we can con- 
clude that the 7-velocities are a side-effect of the line asym- 
metries and related to the problem of determining line po- 
sitions of asymmetric lines. This view is also supported by 
the high conformity of the slopes qq of the V^A^ relation for 
all sample stars (see Tab. [ij which (almost) agree within 
la uncertainties. 

Although this method of arbitrarily modifying the line 
asymmetry is a rather ad hoc method, it still allowed us 
to gain some insight into the origin of the 7-velocities. In 
particular, it shows (1) that the residual 7-velocities are 
related to the shapes of the spectral hues and, consequently, 
to an intrinsic property of Cepheids and (2) that the main 
physical question is not to understand the 7-velocities, but 
to understand the 7-asymmctries. 

Another important point is that this interpretation pro- 
vides a physical meaning to the 7-velocities and, hence, a 
physical reference. The relation between 7-asymmetry and 
7- velocity allows us to compute the contribution of the dy- 
namical Cepheid's atmosphere to the 7-velocity. We can 
thus really provide a 7-velocity corresponding to the space 
motion of the star itself (V^*), independently of the dynam- 
ical structure of its atmosphere. From Tab. [1] we find an 
average systematic red-shift correction of &o = 1.8 ± 0.2 
km s~^ (averaged over 8 stars) between our physically cal- 
ibrated 7-velocities and the ones found in the GCD, com- 
monly used by the community. Consequently, the K-term 
(blueshifted) found in the literature is not due to the kine- 
matic structure of the Galaxy, but to a bias in the previ- 
ous methods of deriving the 7-velocities (cross-correlation, 
Gaussian fit of the spectral line), most likely due to these 
AyVy linear relations. By using only one metallic line to de- 
rive the 7-velocity, one can make, for instance, an error (or 
find inconsistencies) ranging from -2 to km s~^. Finally, 
our results are thus consistent with an axisymmetric rota- 
tion model of the Milky Way. 



4.2. "f -asymmetries 

Different aspects of the 7-asymmetries must be pointed out. 

— The 7-asymmetries arc different from one spectral line 
to the other. This was the motivation for our Ay V-y re- 
lations. Moreover, Ay can be positive or negative. 

— Given the precision of our data, we find no partic- 
ular link between the line depth (i.e. where the line 
forms within the atmosphere, see Paper II) and the 7- 
asymmetries. Moreover, no relation was found between 
the 7-asymmetries and the wavelength. 

— In order to investigate phase shifts, which might be re- 
lated to the 7-asymmctries, we integrate radial veloc- 
ities over time (after correcting by the corresponding 
Vy(i)) and obtain the radius as a function of the pulsa- 
tion phase for different spectral lines. Such phase shifts 
are observed clearly only for R TrA and RS Pup (see 
Fig.SD. 

— One can see in Fig. that the difference in asymme- 
try between the Fe I 6024. SA (dot-dashed line) and Fe I 
4896. 4A (dashed line) asymmetry curves is surprisingly 
almost constant with the pulsation phase. Even it is 
not a general rule (see solid line of Fig. [2b), this effect 
is frequent in our data, and seems to be important to 
understand the meaning of the 7-asymmetries. 

— Although the A{(j)) and RVc{4>) curves vary in phase, 
there is no strict correlation between the line asymme- 
try A{(j)) and the radial velocity i?X4(0) as a function of 
the pulsation phase (j) (see Fig.[5J:). The slope and shape 
of the RVc~A curves depend on the limb-darkening, the 
spectral line width, and the rotation velocity. However, 
the loop and the shift in asymmetry (7-asymmetry) are 
mainly due to the dynamical structure of the Cepheid 
atmosphere. They cannot be explained by static geo- 
metrical models of pulsating stars (see Paper I for de- 
tails). 

— In the stellar rest frame (defined by V^+), one has for a 
given spectral line and at a specific pulsation phase (/), 
RVc{<l>) = while A{(f)) ^ or, on the contrary, 
Ai(t>) = but RVc{(f>) ^ 0; i.e. i?Vc(0) does not 
imply A{(f>) ~ 0. 

— Interestingly, we find a relation between the 7- 
asymmetries and the pulsation period of Cepheids. In 
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particular, if wc define {A^) as the average of the 
quantity over all spectral lines for a given star, then we 
find the following linear relation 

{A-y) ^ [-1.6 ±0.3] logP+ [6.7±0.9], 

where P is the period in days of the Cepheid. The re- 
duced we obtain is 6. This relation (hereafter (A^)) 
is represented by Fig. [5] We already presented such a 
relation in Paper I for the Fel 6056. 005A spectral line. 
This new relation can be considered as a generalization 
over all spectral lines of our sample. 

What could explain the observed 7-asynimetries? To an- 
swer this question, we investigate three different hypothe- 
ses: (1) the limb-darkening variation with time and within 
the spectral line, (2) velocity gradients, and (3) the relative 
motion of the line-forming region compared to the corre- 
sponding mass elements. 

(1) The time- and wavelength dependence of the limb- 
darkening within the spectral lines does not explain why 
the line asymmetry is not zero when the radial velocity is 
zero. When for instance the star is at minimum or maxi- 
mum expansion (i.e. radial velocity is zero), whatever the 
limb-darkening distribution might be, if there are no dy- 
namical effects within the atmosphere, then the spectral 
line should be symmetric. This is clearly not what we ob- 
serve. Consequently, the time- and wavelength dependence 
of the limb-darkening may have some effects on the varia- 
tion of line asymmetry over the pulsation phase, but it can 
not be responsible for all properties mentioned above. 

(2) Velocity gradients in the atmosphere can indeed 
cause an asymmetric spectral line. Indeed, line asymme- 
try A{(j)) and radial velocity RVc{(j)) basically vary in phase, 
even if there is not a strict correlation. This is thus a strong 
indication that radial velocity gradients are, to some degree, 
responsible for the line asymmetries in these stars. However, 
they cannot explain, for instance, the systematic shift in 
asymmetry (almost constant with the pulsation phase) ob- 
served between the 6024. SA and Fel 4896. 4A spectral hue 
(Fig. [2). Butler et al. (1996) introduced a number of veloc- 
ity gradients in their model. They found that "the effect 
of this velocity gradient is to reduce the amplitude of the 
pulsational velocity curve at optical depth of r = | by 20% 
and to decrease the 7-velocity by 2 km s~^relative to the 
standard assumption of a comoving atmosphere" . However, 
the main differences they obtain in their velocity curves are 
near extrema. For instance, in their Fig. 11, all curves (with 
and without velocity gradient) vanish at the same pulsation 
phase. This is not what we observe in our Fig. [2^, where 
systematic shifts are obtained. In their study, Butler et al. 
(1996) applied a "closure" constraint: "It is assumed that 
each layer of the stellar atmosphere returns to its starting 
position after a pulsation cycle." Thus, let us now discuss 
the last hypothesis: 

(3) The relative motion of the line-forming region with 
respect to the corresponding mass elements. When the line- 
forming region moves relative to the background atmo- 
spherical structure, it will also move with respect to the 
background velocity field. Thus, the line experiences an ap- 
parent change in velocity, which has a comparable effect as 
a velocity gradient, but with a lower intensity. This explains 
why line-forming regions do not comply with the path con- 
servation, as discussed in the introduction. However, this 



cannot explain, for instance, the systematic shift in asym- 
metry (or in velocity) observed between the 6024. SA and 
Fe I 4896. 4 A spectral lines. Indeed, it would mean that the 
line-forming regions corresponding to these two spectral 
lines have a systematic relative motion iV^) compared to 
the background, whatever the pulsation cycle considered. 
Either there are cycle-to-cycle differences in the path of 
the line-forming regions (which we cannot confirm with our 
data) or one should invoke another physical explanation. 

5. Conclusions 

We found 7-asymmetries varying from one spectral line to 
another, as well as a global dependency on the period of 
the star or, correspondingly, on the spatial extension of the 
Cepheids' atmosphere. Right now, we have no clear physi- 
cal explanations for this effect. Most likely, it results from a 
combination of several effects in the dynamical structure of 
the Cepheids' atmosphere, such as phase- and wavelength- 
dependence of the intensity distribution within the differ- 
ent spectral lines, velocity gradients, non-linear pulsational 
effects, shock fronts, and relative motions between line- 
forming regions (specific to each spectral line considered) 
and the material. 

In order to further investigate the line asymmetries, 
improved numerical models are required. Since Cepheids' 
atmospheres are not in a hydrostatic state but charac- 
terized by pulsational dynamics, one has to perform non- 
linear, time-dependent simulations of the underlying pul- 
sation. Snapshots from this temporal evolution (including 
the velocity field) can then be used to compute a detailed 
frequency-dependent radiative transport. In order to re- 
solve narrow features such as shock fronts or sharp ion- 
ization regions, high spatial resolution, especially in the 
line-forming regions, is needed. Although convective trans- 
port plays only a minor role in the stellar structure of 
Cepheids in that temperature range (~ 55001^, see Tab. 2 
in Paper II), there can still be considerable convective ve- 
locities (some km s~^). As a consequence, the consistent 
inclusion of the convective velocity field - and of the in- 
teraction of convection with pulsation - in the numerical 
models might be crucial for the computation of line asym- 
metries. Dynamics in the cireumstellar envelope (Kervella 
et al. 2006; and Merand et al. 2006, 2007) could also be 
of importance. Confronting such models with observations 
(spectral line profiles, spatial- and spectral- visibility curves 
from interferometry) may finally lead the way to a complete 
picture of the relevant effects in Cepheids' atmospheres. 

However, even now, we already have a clear evidence 
from our observed linear A^V^ relation that the residual 
7-veloeities (or K-term) seen in Cepheids are the result of 
the dynamical structure of their atmosphere. This provides 
a physical meaning to the 7-veloeities and a physical ref- 
erence: the j-velocity should be zero when the ^-asymmetry 
is zero. This definition of the 7-velocities could be used for 
kinematical studies of the Galaxy, even though this method 
requires high signal-to-noise and high-resolution spectro- 
scopic observations. Using only one metallic line to derive 
the 7-velocity can inflict errors (or inconsistencies) ranging 
from -2 to 1 km s~^. 
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Fig. 4. Radius (in solar radii) as a function of the pulsation 
phase for spectral hues of our sample. The radial velocity 
curves have been corrected from their 7-velocities in order 
to investigate phase shifts. Important differences are found 
from one line to the other, which are mainly due to the 
velocity gradient within the atmosphere. Actually, there is 
a linear relation between the amplitude of the radius vari- 
ation and the line depth (see for instance Fig. 4 of Paper 
II). This indicates why it is essential to use dynamic pro- 
jection factors. A phase shift is observed only for R TrA 
and RS Pup. 



N. Nardetto and collaborators: High resolution spectroscopy for Cepheids distance determination III. 9 
Table 2. A-^ and for all spectral lines of all stars. 



Line (A) 


R TrA 


S Cru 


Y a 




p \J 


or 




A^[%] 


V^[km s^^] 


A^[%] 


VT,[km s-i] 


A^, [7o] 


V7[km s \ 


A^ [%J 


A r n — ll 

V^^km s J 


Fel 4683.560 


0.40±3.oi 


0.02±o.53 


3. 12^2.41 


— 0.68±0.56 


-4.49±2.58 


0.61±0.43 


-0.44±o.43 


-0.21±o.2i 


Fel 4896.439 


— 2.58±7.39 


0.84±o.55 


-4.80±5.76 


0.80±0.58 


0.82±4.23 


1.58±0.46 


--0.90±0.84 


0.70±0.23 


Fel 5054.643 


13.95±5.75 


— 1.11±0.60 


6.55±4.3o 


-0.86±o.6i 


8.99±4.8i 


0.11±0.47 


4.07±o.8o 


-0.56±o.2o 


Nil 5082.339 


11.19±1.95 


— 1.17±0.66 


7.60±i.74 


— 1. 12^0.61 


8.12±1.77 


— 0.80±0.49 


6.77±o.43 


-0.83±o.i9 


Fel 5367.467 


3.81±0.55 


-1.04±o.78 


6.31±0.58 


— 1.02±o.7S 


1.68±0.49 


-0.74±o.65 


3.82±o.i6 


-0.81±0.27 


Fel 5373.709 


1.78±1.73 


-0.50±0.73 


3.63±i.64 


— 0.69±o.70 


— 3.02±i.54 


0.17±0.66 


3.74±o.38 


-0.68±0.25 


Fel 5383.369 


4.66±o.44 


— 1.25ib0.84 


7.61±0.46 


— 1.42±0.88 


3.90±0.39 


— 1.16±0.66 


6.65±o.i4 


— 1.03±0.28 


Till 5418.751 


5.75±o.72 


-0.35±0.79 


5.75±o.65 


0.02±o.7i 


5.59±o.6i 


-0.24±o.59 


4.85±o.i8 


0.13±0.25 


Fel 5576.089 


2.39±o.7o 


— 0.61±0.74 


4.67±o.68 


-0.94±0.84 


0.77±o.6i 


-0.77±o.6i 


3.45±o.i6 


-0.89±o.30 


rei OoDZ.oOo 


1.96±i.u2 


0.09±0.84 


-0.18±i.ii 


0.12±o.85 


— i.^O±0.96 


U.Ui±0.66 


1 on , , . 

i.ZU±0.25 


— U.U4±o.34 


Fe T 6024 05S 


5.70±0.77 


— 1.11±0.89 


6.98±o.77 


— 1.42±o.90 


4.58±0.67 


-1.26±o.7i 


7.58±o.2i 


-1.45±0.25 


Fel 6027.051 


4.68±i.59 


-0.52±o.78 


6.59±i.49 


-0.99±o.8i 


3.71±i.3i 


-0.56±o.6i 


5.38±o.33 


-1.00±0.27 


Fel 6056.005 


3.02±i.53 


-0.07±o.82 


0.60±1.53 


— 0.22±o.80 


2.33±i.4o 


-0.44±o.65 


-0.43±o.35 


— 0.38±o.30 


Si I 6155.134 


-0.70±1.82 


0.31±o.6i 


1.30±1.67 


-0.16±0.63 


-1.01±2.12 


0.29±0.47 


2.85±o.38 


— 0.53±0.22 


Fel 6252.555 


5.19±o.68 


-0.93±o.93 


4.83±o.7o 


— 0.77±o.84 


5.99±o.6i 


— 1.30±o.70 


4.29±o.i6 


-0.58±0.24 


Fel 6265.134 


3.63±i.io 


-0.63±0.83 


3.75±i.i2 


— 0.80±0.89 


1-11±0.92 


-0.70±0.67 


3.26±o.2i 


-0.80±0.28 


Fel 6336.824 


8.99±i.oi 


-1.45±0.89 


7.59±i.o3 


— 1.46±0.84 


6.62±0.88 


— 1.75±0.65 


6.70±0.23 


-1.44±0.24 


averaqe 


4.34±i.87 




4.23±i.63 




2.61±1.52 




3.69±o.33 




Line ( A j 


C Gem 


RZ Vel 


t Car 


RS Pup 




A^[%] 


V^[km s^^] 


A^[%] 


V^[km s^^] 


A^[%] 




A.,[%] 


VT.[km s^^] 


Pel 4683.560 


0.33±0.46 


-0.44±o.22 


— 6.01±io.oo 


0.24±o.99 


— 1.78±0.49 


— 0.32±o.32 


-5.97±3.i8 


0.64±o.63 


Fel 4896.439 


— 7.95±o.92 


0.96±0.24 


— 5.89±io.86 


0.99±i.o2 


2.10±0.86 


0.33±0.33 


— 3.67±3.32 


1.35±0.75 


Fel 5054.643 


3.30±o.8o 


-0.50±o.22 


3.15±7.02 


— 0.22±i.oi 


2.41±o.78 


-0.56±o.3o 


0.22±2.75 


0.01±o.61 


Nil 5082.339 


6.51±0.44 


-0.62±o.2o 


6.75±2.37 


-l.OOii.oi 


4.42±o.64 


-0.76±0.27 


6.24±i.65 


— 0.53±0.52 


Fe I 5367.467 


3.95±o.i9 


-0.43±o.30 


1.00±0.5S 


-0.09±i.5i 


— 0.87±0.27 


0.13±o.4o 


-0.03±0.39 


0.01±l.l3 


Fel 5373.709 


3.56±o.4i 


-0.59±0.28 


1.51±1.87 


0.03±i.3o 


1.73±0.49 


-0.53±0.38 


-0.04±o.99 


0.15±0.94 


Fel 5383.369 


6.53±o.i6 


-0.84±0.33 


4.01±0.47 


-0.79±i.32 


5.75±o.23 


— 1.06±0.43 


8.66±o.3i 


-0.71±i.o7 


Till 5418.751 


4.75±o.2o 


0.45±0.28 


4.14±o.57 


0.17±1.43 


-0.62±o.3i 


l-07±0.44 


4.65±o.45 


— 0.52±o.96 


Fel 5576.089 


2.39±o.i9 


-0.25±0.35 


-1.49±0.74 


-0.37±i.36 


-0.03±o.2i 


0.04±o.5o 


— 3.67±o.47 


0.30±i.ii 


rei ooDz.oOo 


-0.56±o.3o 


0.35±0.37 


-2.89±i.2i 


0.57±i.46 


— 2.42±o.36 


0.43±0.49 


— 6.09±0.72 


0.73±i.ii 


Fel 6024.058 


6.80±0.24 


— 1.01±0.32 


5.83±o.78 


— 1.23±i.4i 


7.54±o.3i 


— 1.24±o.44 


8.42±o.48 


-1.66±0.97 


Fel 6027.051 


4.19±0.37 


-0.43±0.32 


4.32±i.78 


-0.21±i.32 


2.73±o.42 


— 0.38±0.45 


2.90±i.o5 


— 0.37±0.88 


Fel 6056.005 


-2.40±o.4i 


0.15±o.34 


-6.75±i.79 


0.98±i.58 


-7.00±o.5i 


0.68±o.45 


— 8.53±o.85 


0.53±i.o5 


Si I 6155.134 


1.74±o.43 


-0.09±0.25 


-1.35±2.76 


0.71±i.2i 


0.44±o.64 


0.08±0.36 


-1.49±i.i4 


0.26±0.78 


Fel 6252.555 


3.22±o.i9 


0.07±0.27 


0.81±0.75 


— 0.16±i.4o 


0.77±0.27 


0.22±o.34 


1.66±0.47 


— 0.93±0.84 


Fel 6265.134 


2.20±0.24 


— 0.38±0.35 


-0.13±i.3o 


— 0.21±i.44 


1.14±0.23 


-0.07±0.46 


0.23±o.60 


0.70±i.i3 


Fel 6336.824 


4.70±0.27 


-0.87±0.26 


2.81±i.o9 


-0.88±i.24 


3.22±o.35 


— 0.59±0.35 


3.78±o.77 


— 0.92±0.69 


average 


2.54±o.37 




0.58±2.7O 




1.15±0.43 




0.43±i.i5 






1 10 100 
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Fig. 5. averaged (without weighting) over all lines is 
given as a function of the pulsation period of the star. The 
plot is a generalization of the Fig. 14a of Paper 1. 



